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Hot Super Earths: disrupted young jupiters? 
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ABSTRACT 

Recent Kepler observations revealed an unexpected abundance of "hot" Earth-size to 
Neptune-size planets in the inner 0.02—0.2 AU from their parent stars. We propose that 
these smaller planets are the remnants of massive giant planets that migrated inward 
quicker than they could contract. Wc show that such disruptions naturally occur in 
the framework of the Tidal Downsizing hypothesis for planet formation. We find that 
the characteristic planet-star separation at which such "hot disruptions" occur is i? « 
0.03 — 0.2 AU. This result is independent of the planet's embryo mass but is dependent 
on the accretion rate in the disc. At high accretion rates, M ^ 1O~^M0 yr~^, the 
embryo is unable to contract quickly enough and is disrupted. At late times, when the 
accretion rate drops to M ^ 10^® M0 yr~^, the embryos migrate sufficiently slow to 
not be disrupted. These "late arrivals" may explain the well known population of hot 
jupiters. If type I migration regime is inefficient, then our model predicts a pilc-up of 
planets at R ~ 0.1 AU as the migration rate suddenly switches from the type II to 
type I in that region. 



1 INTRODUCTION 



realised (|Vorobvov fc Bas ul 1200^; iBolev et al.] |2010| ). An- 



Stan dard proto-planetary disc models (| Chiang fc GoldreichI 
1 19971 ) show that the inner ~ 0.1 AU region is too hot to 
allow existence of small solid particles there. Thus planets 
should not be able to grow there. Yet observations of nearby 
solar type stars show that many of them do host planets in 
that inhospitable to planet formation region. The very first 
exoplanet to be convincingly detected had the separation 
of -R 0.05 A U to its p arent st ar and had a mass of about 
that of Jupiter (|Mavor fc Qucloal 19951 ) . Such gas giant plan- 
ets circling their parent stars in a close proximity (7? 0.1 
AU) are now called "hot jupiters". It is believed that they 
are explained by the inw ard radial drift (migration) of the 
planets born further out (|Lin et al.lll996l ). 

The Kepler miss ion has recently pr oduced a number 
of surprising results jBorucki et afllioill ). one of which is 
that there is an even greater number of smaller planets, 
e.g.. Earth-size to Neptune-size, in that region. It is sim- 
ilarly obvious that these smaller "hot" planets also had 
to be brought there from further out by an inward ra- 
dial migration. One practical difficulty in testing this idea, 
though, is that the migration of smaller planets is ex- 
pected to occur in the poorly und erstood "type I" regime 
(|Paardekooper fc Papaloizoul |2008| ) . This is different from 
the better under stood "type H" regime b y which the giant 
planets migrate (|Lin fc Papaloizoul 19861 ). This theoretical 



difficulty leads to a large range in uncertainties in the pre- 
dictions of the de tailed Core Accr etion model calculations 
(e.g., see Fig. 5 in llda fc Linll2008l ). 

Recently, the key importance of the radial migration 
of the earliest gas condensations formed in the massive 
proto-planetary discs by the gravitational instability was 



alytical e stimates (INava kshin 2 010al) and numerical sim- 
ulations IV orobyov fc Basul 120061 . l2010l : iBolev et al.] |2010| : 
ICha fc Navaksh in 2010) show that these condensations can 
migrate all the way from their birth-place in the outer 
R ~ 100 AU disc into the inner ~ few AU disc and be 
disrupted there during the earliest massive disc phase {t ^ 
few X 10^ yrs, typically). 

It was pointed out bv lBolev et al.l (|2010f ) and lNavakshinI 
(|2010al ) that this migration-and-disruption sequence yields 
an unexplored way of forming terrestrial like planets. If dust 
grows and sediments in the centre of the clump and forms 
a solid density core there, then tidal disruption of the gas 
clump may leave a solid core - an Ear th-like proto-planets 
(note the connection to earlier ideas of iMcCrea fc Williamj 
1 19651 : iBosdl 19981 : IBoss et al.ll2002l ) . iNavakshinI (|2010d lbl) used 
a simple spherically symmetric radiation hydrodynamic code 
with the dust grains as a second fluid to delineate the con- 
ditions when such a mechanism for the solid core growth 
can work. Based on the potential promise of these ideas, 
INavakshinI (l2010a^ formulate d the "Tidal Downsizing" (TD 
hereafter; Navaksh"inl l2010al ) hypothesis for planet forma- 
tion. In this picture a partial disruption of a ~ lOA/j gas 
clump (which we also call giant embryos; GEs) leaves a gi- 
ant planet, whereas a complete disruption yields a terrestrial 
like planet. 

In this Letter we continue to assess the potential util- 
ity of the TD hypothesis to planet formation. We note that 
another ing r edien t, muted but not ex plicitly considered by 
iBolev et al.] (|2010l ): INavakshinI (|2010d ]b|]a]). must be included 
in the scheme. To explain it, consider isolated GEs flrst. As 
they contract, their internal temperature increases. At early 
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times the rate of this contraction is controlled by the ra- 
diative cooling rate of the embryo - which is by the rate 
at which the embryo can get rid of the excess energy. How- 
ever, when temperature T2nd ~ 2000 K is reached, molec- 
ular hydrogen disassociates. This process is an efficient en- 
ergy sink, which allows the embryo to contract rapidly - in 
fact collapse hydrodynamically - without the need to ra- 
diate the energy away. The embryo collapse stops only at 
much higher densities, and temperatures as high as lO'' K, 
at which point hydrogen is ionised. The embryo must then 
continue a slower contraction, again regulated by the rate 
at which its energy is radiated away. The collapse is known 
as the "secon d collapse" in the star formation literature 
l|Larson|[l969l). when the "first cores" of masses ~ 5QMj 
collapse ( Masunaea fc Inutsukal I2OO0I I to become "second 
cores", which are the proper proto-stars. 

In the TD hypothesis for planet formation, the second 
collapse may be the last step to making a gas giant planet. 
However, as we show below, this final step is not automat- 
ically successful - planets continuing to migrate rapidly to- 
wards their parent stars may still be disrupted at 7i ~ 0.1 
AU. We suggest this process as a way of fo rming the hot 
Supe r Earths observed by the Kepler mission (jBorucki et al.l 
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Figure 1. The embryo's temperature (in units of 10^ K), density 
(in 10~* g cm""^) and radial size, re, in AU, as a function of 
time, as labelled on the figure. Note the abrupt change near the 
end of the calculation, marking the second collapse, when the 
central temperature reaches T ~ 2300 K. The bump in the central 
temperature near t = 1000 yrs is caused by the formation of a 
~ 5 Mq solid core inside the embryo. 



2 THE SECOND EMBRYO 

In analogy to the star formation literature, we refer to the 
GEs that are mainly molecular, embryo's temperature Te < 
Tbnd , as the "first GEs" ; those where H2 is disassociated are 
termed "second GEs" instead. 



2.1 Contraction and collapse of the first embryo 

To illustrate the main point of this paper, we calculate the 
contraction of a giant embryo with "typical" parameter val- 
ues (e.g., those that appear quite reasonable to u s for a solar 
metalicity disc around a ~ solar mass star; see iNavakshinI 
l2010al l. In particular, the embryo mass is Me = lOMj, the 
normalised dust opacity is fc, — 0.5, and the grain meiss frac- 
tion fg = 0.01. The embryo is in itialised as a first core of 
same mass (see iNavakshinlboiOd ') . 

Figure [1] shows the time evolution of the embryo's cen- 
tral temperature (solid, in units of 10^ K), the gas den- 
sity (dotted, in units of 10^* g cm"'^), and the outer ra- 
dius of the embryo, r^, (dashed, in units of 1 AU). The 
calculation is carried out with an updated versio n of the ID 
:as-du st grains radiative hydrodynamics code of lNavakshinI 
Instead of using an ideal gas equation of state 
with 7 = 5/3, the code now uses the equation of state ap- 
propriate for molecular hydrogen, including disassociation 
and rotational and vibrational degrees of freedom for H2, 
wit h the orto-hydrog en to para-hydrogen ratio fixed at 3:1 
(cf. iBolev et al.ll2006l ) . 

Despite the updated equation of state, the evolution 
of the firs t embryo is quite similar to that of the cases 
studied in INavakshinI (|2010bl ). This may not be entirely 
surprising given the similar insensitivity of the first (gas) 
cores t o the equation of state as fou nd bv iMasunaga et al.l 
(|l998l ): lMasunaga fc Inutsukal (|2000l ). The embryo contracts 
and heats up, whereas dust grains grow. By time t ~ 1000 



yrs, the grains increase in size to about 20 cm. Their den- 
sity exceeds that of the gas in the centre of the embryo; 
they become self-gravitating and form a solid core of mass 
Mc ~ 5M0. In Figure [Jl the solid core formation is notable 
by the bump in the central temperature. After the core for- 
mation, the central region becomes hotter than grain vapor- 
isation temperature of ~ 1400 K, evaporating the grains, 
and th us terminating further core growth (see INavakshinI 
l2010bl . for details on this negative feedback loop). The cen- 
tral region also expands slightly. Most of the GE is however 
unaffected by the solid core in this case, and the curves re- 
sume their otherwise monotonic behaviour a few hundred 
years later. 

At t ~ 6.5 X 10^ yrs, the GE goes through the second 
collapse when the central temperature exceeds about 2300 
K. The embryo radius drops rapidly, while the density and 
the temperature increase strongly. The first embryo becomes 
the second in our terminology. 



2.2 The contraction of the second embryo 

Our radiation hy drodynamics code is not well suited 
(|Navakshinll2010bl ') to follow the long cooling and contrac- 
tion of the second embryo. However, there is a body of 
work on the hydrostatic contr action of very low-m ass stars 
and Jupiter-mass pla nets (e.g.. lGrossman fc Graboskc. 1973) : 
iGraboske et al.lll975l ). which allows us to describe the pro- 
cess with a good degree of confidence. Let t2 be the time 
of the second collapse of the embryo. As is well known, the 
embryo spends the initial contraction stages on the Hayashi 
track. During this phase, the outward energy transfer is 
dominated by convection. The effective temperature, T^g 
is almost constant at logTc ff ~ 3.1 — 3.3 [see Figs 1 in 



I Grossman fc Graboskd (jl973l ): IGraboske et all (jl975l l]. The 

luminosity of the embryo is then 4,-Kr1aT*ff. 
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With this simple cooling model, we solve for the evolu- 
tion of the second GE radius, re{t), 



— — = Anr^aBi^a 
at 



(1) 



where the energy of the GE is defined with the positive sign 
as -Ege ~ GMe/2re. A trivial integration yield^ 



rl 



l + Arl(t-t2) 



(2) 



Here r2 is the GE radius at the time of the second collapse, 
t2, and A = 2^-RaBT*s/{GMl). This model does not take 
into account the electron degeneracy pressure in the GE, but 
that becomes important only after t ^ 10*' yrs, and would 
make disruption of the planets even more likely. 

Evaluation of equation ((2| shows that at times {t — 
t2) ^ 100 yrs, the initial radius of the embryo is quickly for- 
gotten, and the contraction proceeds as 



1/3 



^A{t-t2)^ 

The second embryo's density as a function of time is 



3A'hA{t - 12) 
47r 



■{t-t2) 



(3) 



(4) 



In terms of absolute values, the initial density of the second 
embryo is at least 10~^ g cm~^ (cf. Fig. 1) and rises with 
time according to 



«^in-4 -3 {t-t2) lOMj 
p.W = 6xlO gem 

where we set logTefi = 3.1 for certainty. 



(5) 



~ ^ - yrs 



M 



(7) 



Now, this time scale is to be compared with the Kelvin- 
Helmholtz time scale for isolated giant planets, which is of 
th e order of tKH ~ (10 "* — 10^) yrs (e.g., see Figs. 2 & 3 
of iGraboske et al.|[l975l ). In fact, any addit ional effects we 
can think off, e.g., stellar irradiation (e.g., ICameron et al.l 
Il982h . tidal heating, energy stored in the rotation of the 
GEs, accretion of planetesimals, etc., should only increase 

We can now estimate the embryo-star separation at 
which the disruption occurs, 7?hot, e.g., where pe = p{Rhot- 
For this we note that {t — t2) in equation Q should be of 
the order of imig. Indeed, this is the time it takes the em- 
bryo to migrate inwards, and thus the embryo's age should 
be comparable to that. Entering tmig instead of {t — 12) in 
equation Q, we find that the dependence on the embryo's 
mass cancels out, and we arrive at 



-Rhot — 



GM,M 
367raT4j 



1/3 



= 0.12 AU 



M 



10-5 Meyr- 



1/3 



,(8) 



where we set logToft =3.1 and retained the dependence 
on M only. We term this inner disruption "hot" to dis- 
tinguish from the more distant cold disruptions. Equation 
[8] shows that a second GE entering the inner ~ 0.1 AU 
in an early disc phase, when M is very high, is likely to 
be tidally disrupted. On the other hand, "late arrivals", 
when M ^ 10~* M0 yr~^, may survive and contract into 
hot jupiters as the disc runs out of mass, p resumably due to 
photo-evaporation jAlexander et al.ll2006f ). 



2.3 Vulnerability of young massive planets 

A distance R away from the star, the tidal density is 



Pt 



27ri?3 



S cm- ^i?- , 



(6) 



where M* is the stellar mass, and Rau ~ R/^ AU. A 
giant embryo that migrated to a distance R is disrupted 
if the embryo's density pe ^ p%- The first and the second 
embryos have considerably different internal densities. The 
"first embryos" ar e characteristically disrupted at _R ^ 2 AU 
(|Navakshinl2010al ). We term these disruptions "cold" and do 
not consider any further here. 

To evaluate the second GE's fate, we need to estimate 
its migration rate. We assume that there is only one such 
GE in the inner R few AU at any one time. Due to its 
significant m ass, the GE migrates in the type II regime. 
As shown bv llvanov et al.l (| 19991 ). the migration time, imig, 
is shorter than but is comparable to the "accretion time", 
fa = Me/M, where M is the accretion rate in the disc. In 
the TD model, the planets are born in the very early gas- 
rich phase of the disc when the proto-star is far from its final 
mass. The accretion rate during that phase is M ~ M^/tfi, 
where ~ 10^ yrs is the free-fall tim e for typical in terstellar 
molecular gas cores of ~ 1 M© mass (|Larsonlll969l '). Thus, 



^ A similar model can be used to descr ibe the contrac tion of 
young Mt ~ Mq proto-stars as well (e.g., [Cameron|[l995l V 



3 AN ILLUSTRATIVE MODEL 

We now present two approximate example calculations that 
are complete in the sense that they start off with the em- 
bryo's birth at i? = 100 AU and they end with the embryo 
arriving in the innermost disc. We use the model embryo 
calculated in jj2.1 and 2.2 , and an approximate radial migra- 
tion model of iNavakshinI (|2010al 'l. We consider two opposite 
limiting cases to illustrate the points made in 12.31 



3.1 A disrupted hot jupiter 

In the first case, presented in Figure [51 the accretion disc is 
massive and the accretion rate is high. In particular, the ini- 
tial mass of the star is M, = 0.5 M©, the doubling timescale 
for the star, tdb = 10^ yrs. This yields an accretion rate 
through the disc of 5 x 10"^ M© yr^^. At t = fdb, the as- 
sembly of the star is assumed complete (its mass reaches 1 
Mq), and the disc torques are abruptly removed. 

For simplicity we assume that the mass of the GE is 
constant until it is disrupted dynamically. 3D simulations 
of embryos migrating in gas discs show destruction of the 
embryos in a matter of t wo or three orbits (|Bolev et al ]|2010l : 
ICha fc NavakshinlbOldh once they fill their Roche lobes. 

The upper panel of Figure [2] shows the radial location 
of the embryo. The embryo migrates into the inner ~ 0.1 
AU region of the disc in less than 10'' yrs due to the high 
disc accretion rate. The middle panel shows the evolution of 
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the Hills radius, vh, and the embryo's radius, re- is ini- 
tially one to two orders of magnitude smaller than the Hills 
radius, especially right after the second collapse. The lower 
panel shows the embryo's temperature evolution, which is 
comprised of two parts: the first one calculated as in §2.1 
(and shown in Figure 1) until the second collapse occurs, at 
which point we switch to the model of §2.2, with T found 
through the virial relation of the embryo. 

The "hot" disruption of the embryo occurs at time t ~ 
6 X 10'' yrs, at the embryo-star separation of R — 0.06 AU. 
Note that disruption of the gaseous envelope leaves a solid 
core behind since the solid core formation occured much 
earlier. The solid core migrates slightly to R ^ 0.045 AU by 
t = tdb- This last bit of radial m igration is via type I and may 
in reality be far less efficient ijPaardekooper fc Papaloizoul 
l2008l : llda fcTinll2008h . 
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3.2 A tidally stable hot jupiter 

In the other example we take a proto-star that is closer 
to its final mass, A4, = 0.8 Mq, accreting at a rate of 
10"'^ Me yr"' for 2 x 10*^ yrs. Figure |3] shows that in this 
case the embryo migrates much slower. By the time it ar- 
rives in the innermost disc, the embryo radius, re, dropped 
to about 2 X lO"'^ AU. The embryo is too compact to be 
disrupted unless it migrates even further to ~ 0.01 AU (it is 
likely to stall before that due to the magnetospheric cavity 
in the inner disc). 



4 DISCUSSION 

We suggested here that massive gaseous proto-planets in the 
"second" configuration, migrating rapidly inward, are tidally 
disrupted in the inner ~ 0.1 AU. If these protoplanets con- 
tain solid cores, then the cores remain unaffected by the dis- 
ruption as their densities are much higher. For solid cores 
more massive than ~ 20 , parts of the g aseous envelopes 
may be retained as well ()Navakshinll2010b ). Therefore, the 
hot disruption of gas giant planets could in principle result 
in the production of purely solid and also solid plus gas en- 
velope planets similar to those found in the Kepler data. 
Further modelling of the regions closest to the solid cores 
inside the embryos is needed for a more quantitative state- 
ment. 

The outcome of a hot jupiter arriving in the inner 0.1 
AU strongly depends on the migration rate of the latter and 
its age. If the proto-planet is "old", e.g., ^ 1 Myrs, then it 
is too dense to be disrupted near the star. The observed hot 
jupiters in our interpretation are the planets that arrived in 
the inner 0.1 AU somewhat late, when they were already 
compact, and when the disc was running out of material. 

Some of the Kepler Neptune and Earth-size planets 
could have parted with their gaseous envelopes in the cold 
disruptions at a few AU, and then migrated inward in the 
disc. This would make for a second channel of making these 
planets in the inner disc. However, the location of hot dis- 
ruption is quite well defined, R ~ 0.03 — 0.2, due to the 
strong dependence of the tidal density on R. Therefore we 
predict a pile-up of smaller planets at those radii as long as 
the actual type I migration rate of the disruption remnants 
is far smaller than the theoretical type I migration rate (cf. 
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Figure 2. Formation of a "hot Super Earth" starting from a giant 
embryo of 10 Jupiter masses. Upper panel: Radial position of 
the embryo versus time. The GE starts off at R = 100 AU. The 
vertical dashed line marks the time when a solid core forms inside 
the embryo. Middle: Radius of the embryo, re, and the Hills 
radius, rjj. The GE gaseous envelope is assumed to be completely 
disrupted when the two sets of curves meet. Lower: Temperature 
of the first cores as a function of time. Note the strong jump at 
the second collapse. 
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cussions. Richard Alexander is thanked for useful comments 
on the draft of the manuscript. 
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Figure 3. Same as Figure [2] but for the disc accretion rate of 
1O~^M0 yr~^. The embryo is not disrupted in the inner disc 
since it has had enough time to contract to much higher densities. 



Ilda fc Linlliooi ). Such a pile-up may be testable with the 
current exoplanet data. 
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